Systematic study of (p,7) reactions on Ni isotopes 
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A systematic study of the radiative proton capture reaction for all stable nickel isotopes is pre- 
sented. The results were obtained using 2.0 - 6.0 MeV protons from the 11 MV tandem Van de 
Graaff accelerator at the University of Notre Dame. The 7 rays were detected by the NSCL SuN 
detector utilising the 7 summing technique. The results are compared to a compilation of earlier 
measurements and discrepancies between the previous data are resolved. The experimental results 
are also compared to the theoretical predictions obtained using the NON-SMOKER and SMARAGD 
codes. Based on these comparisons an improved set of astrophysical reaction rates is proposed for 
the (p,7) reactions on the stable nickel isotopes as well as for the ^''Ni(p,7)^^Cu reaction. 

PACS numbers: 25.40.Lw, 24.60.Dr, 26.30.Ef 
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I. INTRODUCTION 

During the past several decades great progress has 
been made in understanding nucleosynthesis. Many open 
questions remain nevertheless, such as the production 
mechanism of nuclei on the proton rich side of the val- 
ley of stability. In particular, there are proton-rich iso- 
topes, the so-called p nuclei, that cannot be produced by 
neutron capture processes, as they are shielded from /3 
decay by the valley of stability A number of pro- 
cesses and sites have been suggested for the production 
of these nuclides but each has its own problems Q. The 
currently favored production process, photodisintegra- 
tions in the 0/Ne layers of massive stars during their 
supernova explosion is the so-called 7 process. This pro- 
cess reproduces the bulk of the p nuclei with the excep- 
tion of the lightest ones, with mass numbers A < 100 
d, 01 • Some deficiences have also been found in the re- 
gion 150 < A < 165 Hi]. Different seed abundances, 
not encountered in massive stars, may allow a 7 process 
to also produce light p nuclides. Thermonuclear explo- 
sions of mass accreting White Dwarfs (the single degen- 
erate model for type la supernovae) have been suggested 
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as alternative sites but initially no light p-nucleus pro- 
duction could be achieved |7Hl0l|. Recent simulations, 
however, found that light p nuclei are produced in suffi- 
cient amounts, due to improved hydrodynamic resolution 
of the nuclear burning zones [llj |. 

The original work by [11] suggested to produce proton- 
rich nuclides in a true p process, i.e., by proton cap- 
tures in the H-rich envelope of type II supernovae. This 
was later shown to be unfeasible T3]- Recently, however, 
another process was found in the deepest, still ejected 
layers of a core-collapse supernova, the z/p process (TJ- 
|l6[, which involves rapid proton captures on nuclei at 
and above Ni. It occurs in explosive environments when 
proton-rich matter is ejected under the influence of strong 
neutrino fluxes. When matter in these ejecta expands 
and cools, nuclear statistical equilibrium mainly com- 
prised of "^He, protons and ^^Ni is achieved at temper- 
atures slightly above 4 GK. Rapid proton captures can 
ensue below about 3.5 GK. Within isotonic chains, (p,7)- 
(7,p) equilibrium is established and the nuclei with the 
lowest capture Q values become waitin g po ints where the 
matter flow through (p,7) is halted 117|. Such nuclei 
typically also show long /3 decay lifetimes, significantly 
longer than the expansion timescale. Without the pres- 
ence of neutrinos, the matter fiow would stop already 
at ^^Ge [1^. However, during the explosion the mat- 
ter is exposed to a large neutrino and antineutrino fiux 
[19I [20j . Due to the high proton abundance, antineu- 
trino captures on free protons can produce free neutrons 
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which allow to bypass the waiting points by (n,p) reac- 
tions [H, [21I, [l^l- The path of the vp process is thus 
defined by the balance between the proton captures and 
their inyerse reactions [13, HI] • Its extension to heavier 
nuclei depends on the processing speed which is giyen by 
the (n,p) reaction rates which, in turn, depend on the 
neutron abundance generated by the antineutrinos, and 
the (n,p) reactivities. Accordingly, variation studies have 
found a strong dependence on uncertainties in (n,p) rates 
but smaller dependence on (p,7) rates because the latter 
are in equilibrium most of the processing time [2l| - |23j . 
Nevertheless, there is a change in the final abundances 
when varying the proton captures rates because these fall 
out of equilibrium at the lower end of the 1.5 < T < 3.5 
GK temperature range, within which vp processing oc- 
curs tlSiill-^S]. 

All the rates used in astrophysical vp process simula- 
tions are based on theoretical predictions. It is therefore 
important to find methods based on experimental data to 
constrain these rates. The initi al p art of the t'p-process 
path follows the iV = Z line [ijllli. It is closest to 
stability at the Ni isotopes. The goal of this work was 
to perform a systematic study of the cross sections for 
(p,7) reactions on all stable nickel isotopes, to test the 
models used for the rate predictions and to provide im- 
proved reactivities to be used in the astrophysical simu- 
lations. Earlier measurements exist for the discussed iso- 
topes. However, several cases show significant disagree- 
ments (up to a factor of five difference in cross section) . 
The current measurements also aimed at resolving those 
discrepancies. 

For each of the targets, the same experimental tech- 
nique was applied and the same experimental setup was 
used, to eliminate any uncertainties arising from appli- 
cation of varying detection techniques. All of the results 
were obtained using the 7-summing technique. All the 
data from previous measurements were retrieved from the 
EXFOR database [2^ . The experimental details are pre- 
sented in Sec.ini New and old data are compared to stan- 
dard NON-SMOKER [25!| calculations as weh as to the 
predictions of the new code SMARAGD vO.Q.Os ^ in 
Sec. mil Based on these results, improved reaction rates 
for ^^Ni(p,7)^''Cu are suggested in Sec. IVl 



II. EXPERIMENTAL DETAILS 

The series of experiments described here was per- 
formed using the 11 MV FN tandem Van de Graaff ac- 
celerator at the University of Notre Dame. The targets 
were irradiated with proton beam in the energy range 
2.0 - 6.0 MeV in 0.2 MeV (0.5 MeV in the case of ^'^Ni) 
steps. The measured energy range was chosen to over- 
lap with the previous measurements and when possible 
extended to cover as much of the Gamow window as pos- 
sible. The beam intensity was monitored using a Faraday 
cup. 

The thickness of the nickel targets was measured us- 
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FIG. 1. Sample 7-sum spectra obtained for the (p,7) reaction 
on each of the nickel targets investigated within this work. 



ing the Rutherford backscattering technique (RBS) at 
the Hope College Ion Beam Analysis Laboratory. The 
properties of the targets are listed in Table HI Apart from 
^^Ni, where the enrichment was 40(5)%, all the targets 
were 95(5)% enriched. 

The 7 rays emitted during the irradiation were de- 
tected by the Nal(Tl) segmented summing detector, SuN 
[27} . SuN was developed at the National Superconduct- 
ing Cyclotron Laboratory, Michigan State University for 
capture reaction measurements utilizing the 7-summing 
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TABLE I. Properties of the Ni targets used during the exper- 
iment. 



Isotope 


Enrichment 


Thickness [mg/cm^] 


Q value [MeV] 


MNi 


40(5)% 


0.270(14) 


7.452 


62Ni 


95(5)% 


1.66(20) 


6.122 


s^Ni 


95(5)% 


0.517(67) 


5.866 


eoNi 


95(5)% 


0.676(90) 


4.800 




95(5)% 


0.943(44) 


3.814 



technique. The 7 rays from the decay cascade of the 
entry state populated during the reaction were summed 
within the detector. Thus the final spectra were com- 
prised predominantly of the sum peak at the 7 energy 
= Ec.m. + Q, where i?c.m. is the total kinetic energy 
in the center of mass system and Q is the reaction Q 
value. The data was recorded using the NSCL Digital 
Data Acquisition System (DDAS). 

The data analysis followed the procedure described 
in detail in [2?|. For a given beam energy, the sum 
peak was fitted with a Gaussian combined with a lin- 
ear background. After subtraction of the background 
fit, the sum peak was integrated within the region of 
{E^—3a,EY,+3a), where a is the standard deviation from 
the Gaussian fit, to obtain the total number of events. 
The hit pattern centroid for the events within the same 
range was used to determine the summing efficiency. 

Sample 7 spectra obtained for each of the nickel targets 
are shown in Fig. [T] As can be seen in Fig. [TJ for each 
reaction the sum peak is clearly observed at the higher 
energy range of the spectrum. It can be observed in Figs 
[Ua)-(d), that for a given beam energy, the position of the 
sum peak changes from target to target reflecting differ- 
ent reaction Q values for various isotopes. The lower en- 
ergy range (E < 4.0 MeV) of the spectrum is dominated 
by the room background (not shown here). The main 
background contribution at higher energies comes from 
the cosmic rays. In two cases (^''Ni and ^^Ni), fluorine 
contamination in the target resulted in a peak at 6.92 and 
7.12 MeV from the ^^F(p, (37)^^0 reaction. This contam- 
ination limited the energy range accessible for the mea- 
surements with the ^^Ni target, as for the beam energy 
below 4.4 MeV the sum-peak energy (i?E = Ec.m. + Q) 
would overlap with the contamination lines. Addition- 
ally, in Fig. [Ija) a peak at 7.2 MeV is present that can 
be related to the neutrons from the (p,n) reaction in the 
target. 

During the analysis, contributions from other isotopes 
were also considered. In particular, for the lightest nickel 
isotopes the (p,7) reaction cross section is an order of 
magnitude lower than for heavier isotopes, thus even a 
small concentration of the heavier isotopes could result 
in a non-negligible sum peak. However, no such contri- 
butions were found. 

The (p,7) reaction cross sections of all stable nickel 
isotopes obtained in this work are listed in Table HT] The 
uncertainty of the cross section values includes the sta- 



TABLE II. Cross section for the (p,7) reactions obtained in 
this work for all stable nickel isotopes. 



Ec.n.. [MeV] 


[mb] 


Ec.n,. [MeV] 


[mb] 




''*Ni(p,7) 


''^Cu 




1.96 


0.66(11) 


3.54 


0.12(02) 


2.16 


1.13(18) 


3.73 


0.20(03) 


2.35 


1.39(22) 


3.93 


0.17(03) 


2.55 


1.03(16) 


4.14 


0.26(04) 


2.75 


0.33(05) 


4.33 


0.22(04) 


2.94 


0.15(02) 


4.52 


0.30(05) 


3.14 


0.24(04) 


4.72 


0.22(04) 


3.34 


0.15(02) 


4.92 


0.28(05) 




-Ni(p,7) 


''^Cu 




2.41 


1.08(14) 


3.40 


1.99(25) 


2.90 


1.379(17) 


3.93 


2.73(34) 




''^Ni(p,7) 


''^Cu 




1.95 


0.34(04) 


2.74 


0.92(11) 


2.14 


0.49(06) 


2.84 


0.89(11) 


2.34 


0.64(08) 


2.94 


0.96(12) 


2.54 


0.75(09) 


3.03 


0.83(10) 




''"Ni(p,7) 


''^Cu 




2.73 


0.30(04) 


4.31 


0.97(12) 


2.93 


0.27(03) 


4.51 


1.37(16) 


3.13 


0.25(03) 


4.71 


0.87(11) 


3.33 


0.50(06) 


4.90 


0.91(12) 


3.52 


0.42(05) 


5.10 


2.37(30) 


3.72 


0.65(08) 


5.49 


2.68(34) 


3.92 


0.70(09) 


5.89 


2.32(30) 


4.11 


1.01(13) 






"«Ni(p,7)"^Cu 


4.30 


0.147(19) 


5.09 


0.266(34) 


4.50 


0.161(20) 


5.29 


0.432(56) 


4.70 


0.134(17) 


5.49 


0.482(62) 


4.89 


0.152(19) 


5.88 


0.378(49) 



tistical uncertainty as well as the uncertainty resulting 
from subtraction of the background under the sum peak. 
The uncertainty in the target thickness of 5% was also 
included, as well as the uncertainty of the summing effi- 
ciency obtained using the method described in [27*] (typ- 
ically of the order of 10% relative uncertainty). 



III. THEORETICAL ANALYSIS ~ GENERAL 
CONSIDERATIONS 

The studied reactions proceed via resonant compound 
formation, either through isolated resonances or a high 
number of overlapping ones, and thus depend on parti- 
cle and radiation widths. From the comparisons of the 
experimental results with the theoretical predictions it 
will become evident that each case has to be analyzed 
separately because the dependence of the cross sections 
on the various widths is different for each case. A simple 
comparison of data and predictions without considering 
these dependencies will be misleading and may not be 
able to identify the underlying deficiency trends and thus 
also the impact on the astrophysical reaction rates cor- 
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FIG. 2. Sensitivity of the reaction cross sections of 
''''Ni(p,7)^^Cu to variations of 7 and particle widths. 

rectly. At temperatures relevant for the vp process the 
reaction rates mainly depend on the proton widths. The 
reaction rate sensitivity on the 7 width increases with 
decreasing neutron number of the compound nucleus as 
its proton separation energy S-p is decreasing. Therefore 
the knowledge of the optical proton+nucleus potential, 
required for the calculation of the proton width, is es- 
sential to constrain the rate in all cases. All sensitivities 
and their definition can be found in [2^. Figures [211S] 
show the cross section sensitivities required in the analy- 
ses of the present work. All cross sections are insensitive 
to variations of a- widths in the shown energy range. 

The data are compared to calculations within the sta- 
tistical Hauser-Feshbach model [1^ . This model assumes 
that a large number of resonances is present at the com- 
pound formation energy and that their individual widths 
can be replaced by average widths. Required is a "suf- 
ficiently" high nuclear level density (NLD) at the com- 
pound formation energy = S-p + Ec.m. [30l - [33 |. Ex- 
cept for ^^Cu and ^^Cu, all NLDs are large enough at Ec 
to apply the statistical model. The NLD is clearly too 
low in 58Ni(p,7)^9Q^^ ^^^^ 60Ni(p,7)^iCu seems to be at 
the borderline. Even when the model is not applicable, 
i.e., when resonance structures show up in the excita- 
tion functions, the statistical model should still be able 
to give an average value for the cross sections. If several 
resonances are contributing, this average value may be 
sufficient for the computation of the reaction rate which 
involves an integration over the cross sections. 

Finally, it has to be considered whether stellar rates 
can directly be inferred from the data or whether addi- 
tional theory has to be invoked. This depends on the 
contribution of reactions commencing on thermally pop- 
ulated, excited states of the target nuclei in the stellar 
plasma. A ground state (g.s.) contribution, i.e., the con- 
tribution of the laboratory cross section to the stellar 
one, of Xo ~ 1 is required to ensure that excited states 
do not contribute |33l - [35| . This is fulfilled for all inves- 
tigated reactions, except for ^^Ni(p,7)^^Cu which shows 
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FIG. 3. Sensitivity of the reaction cross sections of 
''^Ni(p,7)''^Cu to variations of 7 and particle widths. 
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FIG. 4. Sensitivity of the reaction cross sections of 
''^Ni(p,7)''2Cu to variations of 7 and particle widths. 
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FIG. 5. Sensitivity of the reaction cross sections of 
''"Ni(p,7)''^Cu to variations of 7 and particle widths. 
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FIG. 6. Sensitivity of the reaction cross sections of 
^*Ni(p,7)^^Cu to variations of 7 and particle widths. 



Xn = 0.52 at 1.5 GK 

In the following, the reactions are discussed one by 
one, starting with ^''Ni(p,7)^^Cu and going towards the 
more neutron-deficient target nuclei. The cross section 
predictions and width variation studies were performed 
with the code SMARAGD, version 0.9.0s [H, H^- Ad- 
ditionally, comparisons to the standard NON-SMOKER 
cross sections and reaction rates [25), i36| are shown. The 
standard optical proton potential used in those calcula- 
tions is from [37|, w hich is a low-energy modification of 
the potential by [39[ and was provided especially for as- 
trophysical applications. The potential of [391] uses the 
Briickner-Hartree-Fock approximation with Reid's hard 
core nucleon-nucleon interaction and adopts a local den- 
sity approximation. Reparameterizations of [soj using 
more recent data are (40| and its Lane-consistent version 
(4l| . They do not, however, adopt special considerations 
for low proton energies. 

Another modification of [37| , using an increased imag- 
inary part, was developed for an improved description 
of low-energy (p, 7) and (p,n) reactions on intermediate 
nuclei [H, |42h43 |. The present data and calculations 
are not sufficient to distinguish between the two versions 
with and without the modified imaginary part. Due to 
the uncertainties in the description of the widths, similar 
results can be obtained using either potential. Therefore 
only calculations using ^3] shown here. 



IV. RESULTS 

A. ^^Ni(p,7)'*5Cu 

The new results agree well with the previous measure- 
ment by [131 . As shown in Fig. [21 the reaction cross 
sections are sensitive to all widths above the (p,n) thresh- 
old (2.495 MeV) while they quickly become only depen- 
dent on the proton width when going to lower ener- 
gies. A comparison with predictions using various op- 
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FIG. 7. Comparison of measured reaction cross sections of 
''''Ni(p,7)''5Cu to predictions: NON-SMOKER ^] (using 
the potential of [33) and SMARAGD with the potentials of 
[13 (Lej), [U (Bauge), and ^ (Bauge2). The data are from 
this work and from [43| (Sevior). 



tical proton+nucleus potentials (and with the standard 
NON-SMOKER prediction) is shown in Fig. ^1 The 
SMARAGD calculation with the potential by ^37| repro- 
duces the data at the lower edge of the energy range 
in which measurements are available but diverges when 
going to higher energy. Our new data allow a detailed 
study of the cause of this divergence because they extend 
to higher energy. 

A rescaling of the proton width cannot lead to an im- 
provement in the predictions as it would affect the cross 
sections at all energies. Since the neutron width only 
impacts the cross sections above the (p,n) threshold and 
the deviations also occur below it, the only possibility is 
a change in the 7 width. Figure [5] shows how a calcula- 
tion with the standard potential by (37j but a 7 width 
divided by 4 reproduces the data well across the full mea- 
sured energy range. 

The 7 width, in turn, has three dependencies: on the 
low-lying discrete excited states of the compound nucleus 
included in the calculation, on the NLD above the last 
included excited state, and on the 7-strength function. 
Thanks to the new data it is possible to distinguish be- 
tween the three possibilities. The major contribution to 
the 7 width comes from 7 transitions leading to states 
about 2-4 MeV below the compound formation energy 
[4^. Since the proton separation energy in ^^Cu is larger 
than these 7 energies, the low- lying, discrete states only 
play a minor role in this case. A variation of the NLD 
yields different changes in the dependence than varying 
the 7 width in total, as shown in Fig. [HI Therefore the de- 
ficiency in the predicted 7 width stems from the strength 
function used, i.e., the GDR description. This is inconse- 
quential for the astrophysical rate, as it is only sensitive 
to the proton width. 

Within the measured energy range, the experimental 
data is well reproduced by the SMARAGD calculations 
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FIG. 8. Comparison of measured reaction cross sections of 
''^Ni(p,7)''^Cu to predictions using a 7 width divided by a 
factor of 4 (O.257), and NLDs divided by 4 (0.25NLD) and 
10 (O.INLD), respectively. The data are from this work and 
from [13] (Sevior). 



with the renormahzed 7 width but unchanged proton 
width. Adopting this prediction also in the astrophys- 
ically relevant energy range of 0.96 - 2.52 MeV [46^] leads 
to a higher rate than the one of the standard NON- 
SMOKER prediction. The new reactivity is given in Ta- 
ble mil The parameters of a fit in the REACLIB format 
are given in Table IIVI The reactivities across a wider 
temperature range are provided as this is necessary for 
using the table directly in nucleosynthesis reaction net- 
works. It also enables users to make their own fit, if they 
so desire. The data only constrain the rates above 5 GK 
but this is of no consequence for reaction networks, as 
above 5 GK all rates are in equilibrium and rate values 
only have to be provided to avoid divisions by zero. 



FIG. 9. Same as Fig. [3 but for ^^Ni(p,7)''^Cu. The data are 
from this work and from [i^l (Tingwell). 
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FIG. 10. Same as Fig.[7]but for "Ni(p,7f ^Cu. The data are 
from this work and from ^ (Krivo77), ;50] (Tingwell88) . 



B. 



2Ni(p,7)«3Cu 



'Ni(p,7)«2Cu 



The cross section sensitivities of ^^Ni(p,7)^'^Cu are 
quite similar to those of ^'*Ni(p,7)^^Cu but the sensitivity 
to the 7 width is more slowly declining towards lower en- 
ergies. Contrary to the situation with ^^Ni(p,7)^^Cu, the 
SMARAGD calculation using the standard proton poten- 
tial by ^jj is in good agreement with the data across the 
measured energy range. The new data obtained within 
this work confirm the previous measurements ^4 81] , as can 
be seen in Fig. |9l At the high energy end, there is a 
slightly larger deviation between theory and data but 
it is only 10 - 20% which is within the expected range 
of accuracy of a global prediction. Again, there is good 
agreement at the lower end of the measured range, where 
the proton width is determining the cross section. 

The stellar reactivity of the standard SMARAGD pre- 
diction is given in Table IIIII and the fit parameters are 
listed in Table HVl 



With decreasing neutron number, the cross sections 
above the (p,n) threshold depend less on the neutron- 
and proton- widths and primarily depend on the 7 width. 
Below the threshold, its importance decreases again and 
the proton width determines the cross sections in the 
astrophysical energy range. This is also the case for ^^Ni. 

Figure [TU] shows the comparison of various predictions 
to the data, similar as in the previous sections. Two data 
sets from [49,^] for ^^Ni(p,7)^^Cu reaction cross section 
were found in the literature, and a significant discrepancy 
between these data exists. Good agreement between the 
current work and the results from [soj was found. The 
SMARAGD calculation u sing the standard potential by 
(371] describes the data of |49| | well. The data of [S^l and 
those obtained within this work can be reproduced by 
reducing the 7 width by a factor of 0.3, see Fig.[TTJ Con- 
trary to the ^'*Ni(p,7)^^Cu case, however, it is impossible 
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TABLE III. Predicted reactivities of (p,7) reactions on all stable nickel isotopes. 



L J 


"*Ni 






»^Ni 






n^Ni 












5«Ni 






0.1 


5.846 


X 


10-'" 


3.991 


X 


10-'" 


3.534 


X 


10-'" 


2.058 


X 


10-'" 


1.240 


X 


10-'" 


0.15 


2.007 


X 


10-^5 


1.372 


X 


10-1^ 


1.215 


X 


10-15 


7.041 


X 


^0-16 


4.213 


X 


10-" 


0.2 


1.359 


X 


10-1' 


9.316 


X 


10-" 


8.252 


X 


10-" 


4.763 


X 


10-" 


2.832 


X 


10-" 


0.3 


4.728 


X 


10"^ 


3.260 


X 


10-'' 


2.891 


X 


10-^ 


1.657 


X 


10-3 


9.744 


X 


10-" 


0.4 


7.984 


X 


10-" 


5.542 


X 


10-' 


4.922 


X 


10-' 


2.791 


X 


10-' 


1.621 


X 


10-' 


0.5 


3.031 


X 


IQ-^ 


2.118 


X 


IQ-^ 


1.883 


X 


10-5 


1.047 


X 


10-5 


5.930 


X 


10-® 


0.6 


4.793 


X 


10-'' 


3.361 


X 


10-* 


2.989 


X 


10-* 


1.614 


X 


10-* 


8.749 


X 


10-5 


0.7 


4.274 


X 


10"3 


2.999 


X 


IQ-^ 


2.665 


X 


IQ-^ 


1.385 


X 


10-3 


7.092 


X 


10-* 


0.8 


2.548 


X 


IQ-' 


1.783 


X 


IQ-' 


1.582 


X 


IQ-' 


7.860 


X 


10-3 


3.789 


X 


10-3 


0.9 


1.128 


X 


10-1 


7.861 


X 


10"' 


6.959 


X 


10-' 


3.294 


X 


10-' 


1.497 


X 


10-' 


1.0 


3.986 


X 


10-1 


2.761 


X 


10-1 


2.436 


X 


10-1 


1.099 


X 


10-1 


4.725 


X 


10-' 


1.5 


2.809 


X 


lOi 


1.887 


X 


IQi 


1.617 


X 


IQi 


5.973 






2.072 






2.0 


3.314 


X 


lO' 


2.199 


X 


lO' 


1.803 


X 


lO' 


5.900 


X 


lOi 


1.780 


X 


IQi 


2.5 


1.653 


X 


10^ 


1.127 


X 


10^ 


8.745 


X 


lO' 


2.689 


X 


lO' 


7.363 


X 


IQi 


3.0 


5.005 


X 


10^ 


3.669 


X 


10^ 


2.660 


X 


10^ 


8.052 


X 


lO' 


2.051 


X 


lO' 


3.5 


1.100 


X 


10* 


9.013 


X 


10^ 


6.025 


X 


10^ 


1.868 


X 


10^ 


4.487 


X 


lO' 


4.0 


1.941 


X 


10* 


1.830 


X 


10* 


1.115 


X 


10* 


3.666 


X 


103 


8.355 


X 


10' 


4.5 


2.925 


X 


10* 


3.233 


X 


10* 


1.782 


X 


10* 


6.388 


X 


103 


1.385 


X 


103 


5.0 


3.915 


X 


10* 


5.120 


X 


10* 


2.550 


X 


10* 


1.017 


X 


10* 


2.100 


X 


103 


6.0 


5.410 


X 


10* 


9.952 


X 


10* 


4.099 


X 


10* 


2.077 


X 


10* 


3.924 


X 


103 


7.0 


5.814 


X 


10* 


1.486 


X 


10^ 


5.241 


X 


10* 


3.337 


X 


10* 


5.913 


X 


103 


8.0 


5.195 


X 


10* 


1.811 


X 


10^ 


5.717 


X 


10* 


4.359 


X 


10* 


7.547 


X 


103 


9.0 


4.052 


X 


10* 


1.876 


X 


10^ 


5.565 


X 


10* 


4.798 


X 


10* 


8.468 


X 


103 


10.0 


2.895 


X 


10* 


1.707 


X 


10^ 


4.987 


X 


10* 


4.623 


X 


10* 


8.596 


X 


103 



TABLE IV. Fit parameters (in REACLIB format [3^) for the reactivity of (p,7) and its reverse reaction for all stable nickel 
isotopes. The reverse value has to be multiplied by the ratio of the partition functions to obtain the (7,p) reactivity. 





fflo 


ffli 


12 


0.3 


04 


as 




ffle 


^Ni 




7.414078 X IQi 


-1.175453 


2.162346 X IQi 


-1.003880 X 10' 


5.252775 


-3.743412 X 10" 


-i 


4.653832 x lOi 


rev 


9.643699 x IQi 


-8.766878 x Ifli 


2.162346 X IQi 


-1.003880 X 10' 


5.252775 


-3.743412 X 10- 


-1 


4.803832 X Ifli 


•"Ni 




1.173687 X 10' 


-2.367050 


8.350530 X IQi 


-2.120871 X 10' 


1.315014 X IQi 


-8.584441 X 10- 


■ i 


9.604393 X IQi 


rev 


1.396641 X 10' 


-7.341479 X IQi 


8.350530 X IQi 


-2.120871 X 10' 


1.315014 X IQi 


-8.584441 X 10- 


■ 1 


9.754393 x IQi 


^iNi 




9.846000 X IQi 


-1.722409 


5.297264 x IQi 


-1.598366 X 10' 


9.315694 


-6.041871 X 10- 


-i 


7.246493 x IQi 


rev 


1.217359 X 10' 


-6.979359 x lOi 


5.297264 x IQi 


-1.598366 X 10' 


9.315694 


-6.041871 X 10- 


-1 


7.396493 x IQi 


«"Ni 




1.306414 X 10' 


-1.868247 


7.166889 X IQi 


-2.159482 X 10' 


1.421891 X IQi 


-9.256727 x 10" 


-i 


9.242814 X IQi 


rev 


1.529360 X 10' 


-5.757600 X IQi 


7.166889 X IQi 


-2.159482 X 10' 


1.421891 X IQi 


-9.256727 x IQ- 


-1 


9.392814 X lOi 


5«Ni 




1.025240 X 10' 


0.000000 


-1.002597 X IQi 


-1.027998 X 10' 


7.767831 


-5.337708 X 10" 


-i 


3.462199 X IQi 


rev 


1.248178 X 10' 


-3.967033 X IQi 


-1.002597 X IQi 


-1.027998 X 10' 


7.767831 


-5.337708 X 10- 


-1 


3.612199 X lOi 



to decide the reason for the reduction. As shown in Fig. 
[TT| the variation in the NLD leads to indistinguishable 
results from the variation of the full 7 width. A simple 
rescaling of the 7 strength would lead to the same result. 

The comparison of the SMARAGD prediction is com- 
patible with an unchanged proton width and therefore 
we can assume that the astrophysical reactivity is also 
predicted well. It is given in Table|lll]and the REACLIB 
fit coefficients are in Table ITVl 



D. ''"Ni(p,7)«iCu 

The proton separation energy in ^^Cu is only 
4.801 MeV whereas the neutron separation energy is 
11.711 MeV. The data compared to predictions in Fig. 
[12] do not reach the (p,n) threshold and also only barely 
reach the region where the proton width is dominating 
at low energy. The reaction rate would be equally sen- 
sitive to proton and 7 width at the high temperature of 
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this work 1 — ♦ — 1 
Tingwell88 ' — ' — ' 
Krivo77 i-^^ 

SMARAGD(Lej) 

SMARAGD(0.3y) 

SMARAGD(0.3NLD) 



1 1.5 2 2.5 3 3.5 4 4.5 5 



Ec.m. [MeV] 

FIG. 11. Comparison of measured reaction cross sections of 
''^Ni(p,7)^^Cu to predictions using a 7 width multiplied by 
a factor of 0.3 (O.37) and NLDs multiplied by the same fac- 
tor (0.3NLDV The data are from this work and from [HI 
(Krivo77), [50] (Tingwell88). 



10" 




1 1.5 2 2.5 3 3.5 4 4.5 5 



Ec.m. [MeV] 

FIG. 12. Same as Fig. mbut for ^°Ni(p^)^iCu. Data are 
from this work and from [49| (Krivo77), [4g| (Tingwell). 



3.5 GK but rapidly becomes sensitive to only the pro- 
ton width towards lower temperature, including the vp 
process freeze-out temperature of 1.5 GK. 

In the case of the ^°Ni(p,7)^^Cu reaction, the two data 
sets d^, found in literature do not agree with each 
other. The results of this work, as can be seen in Fig. [121 
are in a good agreement with those of [il] and extend to 
higher energies, allowing to better constrain the energy 
dependence. The SMARAGD calculation is close to the 
new data at the upper end of the measured range but is 
closer to the result of [i^ otherwise. 

It proved difficult to obtain the energy dependence of 
the combined data of our measurement and the one by 
(isj . It is only possible by using a scaled NLD without in- 
cluding experimentally known excited states in the com- 
pound nucleus. Figure [T3| shows the theoretical results 
using a NLD reduced by a factor of 5, with and without 




1 1.5 2 2.5 3 3.5 4 4.5 5 



Ec.m. [MeV] 

FIG. 13. Comparison of measured reaction cross sections of 
''°Ni(p,7)''^Cu to predictions using NLDs divided by 5, with 
(0.2NLD) and without (0.2NLD,noexp) excited states for the 
compound nucleus included in the calculation. Data are from 
this work and from ^ (Krivo77), [13] (Tingwell). 



including discrete excited states. With discrete excited 
states, the renormalization of the NLD only affects the 
7 width at higher energy, thus leading to a better re- 
production of the data at higher energy but leaving the 
discrepancy at lower energy. A slightly better energy de- 
pendence can be obtained when applying a factor of 0.3 
to the 7 width and leaving the NLD unchanged. 

This raises the question of the completeness of the ex- 
perimental level scheme. Discrete excited states were 
taken from the 2010 versions of [5l|, [l^l, up to an en- 
ergy of 3.943 MeV. It is a well-known problem in Hauser- 
Feshbach calculations that it is advantageous, on one 
hand, to include low-lying excited states but, on the other 
hand, the included level information has to be complete 
to guarantee accurate predictions. It is often hard to de- 
cide at which excitation energy a cutoff should be made, 
especially in global calculations, in which the cutoff has 
to be implemented through some automated algorithm. 
The difficulty reproducing the energy dependence of the 
measured cross sections encountered here points to such a 
problem. For another recent case and further discussion, 
see [53| . 

The data do not extend very much into the energy 
region where the proton width dominates the energy de- 
pendence. Similar to the other proton capture reactions 
discussed above, it seems that the proton width is well 
described by the optical potential of [37]. Therefore, the 
reactivities and their fit coefficients, shown in Tables IIIII 
and llVl respectively, are based on the SMARAGD calcu- 
lation using this proton potential and the renormalized 7 
width without inclusion of discrete excited states (shown 
as solid line in Fig. [T5)) . 
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FIG. 14. Comparison of the reaction cross sections of 
^*Ni(p,7)^®Cu to predictions. Data are from this work and 
[iil (Krivo77), >4] (Tingwell85), [H (Cheng). 



10"- 



10" 



10- 



10'' 



10'' 




this work 
Tingweil85 
Krivo77 
Cheng 
SMARAGD 
SMARAGD(0.3y) 
SMARAGD(0.3y,Iow exc. stat.) 



3 



[MeV] 



FIG. 15. Comparison of measured reaction cross sections 
of ^*Ni(p,7)^^Cu to predictions using excited states in the 
compound nucleus up to above the proton separation energy 
(SMARAGD), the same but with the 7 width reduced by a 
factor of 0.3 (0.87), and a calculation using a limited set of 
excited states and a reduced 7 width (0.87, low exc. stat.). See 
text for details. Data are from this work and [49l ] (Krivo77), 
[H (Tingwell85), [H (Cheng). 



the fact that SMARAGD is using a more recent version 
of [51I, [ill. Both codes yield almost the same results, 
with a very similar energy dependence which is different 
from the experimentally found one. 

Similar to the ^"Ni case discussed above, better 
agreement between calculation and experiment can be 
achieved by including fewer excited states in the calcu- 
lation. Figure [15] shows how a combination of a 7 width 
reduced by a factor of 0.3 and a limitation of the level set 
to excitation energies i?exc < 2.391 MeV (the J'^ = 9~ 
level) reproduces the mean data over a wide range of 
energies. Considering only the g.s. and neglecting all 
excited states gives the same result, indicating that the 
NLD is predicted well in this confined energy range. Us- 
ing the same 7 width renormalization with the full level 
scheme results in a too low cross section at higher ener- 
gies, indicating the incompleteness of the level scheme. 

The data by [i^ [s^, [s^ and from this work agree 
with each other well although the cross section varia- 
tion structures seem to be different. It has to be real- 
ized, however, that no high-resolution experiment was 
performed, able to resolve the (partially overlapping) in- 
dividual resonances expected for this reaction. Depend- 
ing on the exact beam energy and beam profile, differ- 
ent parts of the same resonances are sampled, leading to 
seemingly different patterns. Therefore, the larger scat- 
ter seen in the data by ^55j] may well be compatible with 
the other data. As mentioned in Sec. IIIIl a statistical 
Hauser-Feshbach calculation can only aim to reproduce 
the average cross sections, which seems to be success- 
fully done by the renormalized calculation shown in Fig. 
1151 For the same reason, however, it may be misleading 
to use the experimental data points to compute a reac- 
tion rate. They do not fully resolve the resonances and a 
simple interpolation will lead to incorrect results. Using 
the theoretical, already averaged values will yield a more 
realistic rate, when lacking further knowledge of the reso- 
nance properties. Therefore, the reactivities provided in 
Table [mi are based on the calculation using the renormal- 
ized 7 width and the limited set of excited states (shown 
as soHd Une in Fig. [15]). The REACLIB fit coefficients 
for these reactivities are given in Table IIVI 



E. 5*Ni(p,7)59Cu 

The NLD at compound formation energy is clearly too 
low for this reaction to expect that all features of the 
cross sections can be described by an average over res- 
onances. Moreover, we encounter a similar problem as 
for ^''Ni(p,7)^-'^Cu regarding the completeness of the in- 
cluded discrete level scheme of ^^Cu. The proton separa- 
tion energy in '^^Cu is Sp — 3.419 MeV. The calculations 
of NON-SMOKER and SMARAGD, shown in Fig. [Mj 
made use of 19 levels (i^^ax = 3.1141 MeV (H) and of 29 
levels (i?max = 4.307 MeV) above the g.s., respectively. 
The levels included below Sp are very similar, despite of 



V. CONCLUSIONS FOR 5''Ni(p,7)"Cu 

One has to be careful in drawing conclusions for the 
stellar ^^Ni(p,7)^^Cu rate from the above trends. It is not 
clear to what extent the conclusions can be applied to the 
doubly magic nucleus with a very low proton separation 
energy of 0.695 MeV. Its reactivity depends on both the 
proton and 7 widths down to 1 GK, as shown in Fig. 1161 
Due to the low Q value, however, it will equilibrate with 
its reverse rate already at lower temperature than most 
other proton captures in the ^-p-process path. 

In the discussion above it was found that the proton 
widths were predicted well in all cases. A renormalization 
of the 7 widths by a factor of 0.3 (0.25 for 64Ni(p,7)65Cu) 
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FIG. 16. Sensitivity of the stellar reactivity of ^®Ni(p,7)"Cu 
to variations of 7 and particle widths, shown as function of 
plasma temperature T. 
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FIG. 17. Ratio of the newly derived rate for ^*'Ni(p,7)'^'^Cu 
and the standard NON-SMOKER rate [l^, shown as function 
of plasma temperature T. 



TABLE V. Predicted reactivity of '^'^Ni(p,7)^'^Cu as function 
of plasma temperature T. 



T [GK] 


A^A (av) 

[cm'^ s~^ mole~^] 


0.10 


9.665 X 10"^' 


0.15 


3.243 X lO"^*^ 


0.20 


2.154 X 10'" 


0.30 


6.990 X 10"" 


0.40 


9.798 X 10"* 


0.50 


2.758 X lO"*^ 


0.60 


3.058 X 10"^ 


0.70 


1.895 X 10"* 


0.80 


7.983 X 10"* 


0.90 


2.571 X 10"^ 


1.00 


6.811 X 10"^ 


1.50 


1.669 X 10"^ 


2.00 


1.032 


2.50 


3.445 


3.00 


8.269 


3.50 


1.628 X 10^ 


4.00 


2.821 X 10^ 


4.50 


4.473 X 10^ 


5.00 


6.653 X 10^ 


6.00 


1.286 X 10^ 


7.00 


2.179 X 10^ 


8.00 


3.319 X 10^ 


9.00 


4.543 X 10^ 


10.00 


5.533 X 10^ 



was required in all cases, except for ^^Ni(p,7)^^Cu. The 
data for the latter reaction was reproduced well by the 
SMARAGD prediction but slightly increasing deviations 
at higher energy may indicate that a suppression of the 7 
width may be required, too. Problems with the NLD only 
arose from incomplete level schemes of low-lying states, 
not from theoretical NLDs. 

Applying the above, an educated guess for an improved 



^^Ni(p,7)^''Cu reactivity can come from a SMARAGD 
calculation using the proton potential by [131, a renor- 
malized 7 width, and no discrete excited states. It has 
to be emphasized once again that such a calculation can 
only yield an average value for the cross sections which 
will only be a good approximation for obtaining the reac- 
tivity if, folded with a Maxwell-Boltzmann distribution, 
they yield a similar value as the actual resonance con- 
tributions. This will not be the case if only few, widely 
separated resonances contribute to the reaction rate in- 
tegral [ssf. 

Table fVl gives the reactivity values obtained with the 
above procedure for the ^^Ni(p,7)^''Cu reaction and the 
REACLIB fit parameters are provided in Table IVTl The 
ratio of the new rate to the previously widely used stan- 
dard NON-SMOKER rate is shown in Fig.HTl In the rel- 
evant temperature range, the new rate is lower than the 
previous standard rate by a factor of about 0.37. How- 
ever, the new reaction rate has a negligible impact on the 
final abundances predicted by the i^p process calculations 
due to the aforementioned equilibrium between the (p,7) 
and (7,p) reactions. 

VI. SUMMARY AND CONCLUSIONS 

In summary, the cross sections of the (p,7) reaction 
for all stable nickel isotopes have been measured. The 
results were compared with previous data found in lit- 
erature. For 64,62,58j^j good agreement with previous 
results was found. For ^^'^"Ni, the new results confirmed 
those of [3 [H3| and disagree with the values of 49]. 
All the results were compared with the standard NON- 
SMOKER calculations and with new predictions of the 
SMARAGD code. New reaction rates were proposed for 
all the reactions studied in the present work. In addi- 
tion, a new estimate of the ^^Ni(p,7)^^Cu rate was de- 
rived based on the comparisons of predictions and data 
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TABLE VI. Fit parameters (in REACLIB format ;36i) for the reactivity of ^*'Ni(p,7)"Cu and its reverse reaction. The reverse 
value has to be multiphed by the ratio of the partition functions to obtain the (7,p) reactivity. 





ao 


ai 


Cl2 


as 


04 


as 


ac, 




7.278475 x 10" 


0.000000 


-4.574355 x 10" 


-3.572965 x 10" 


3.975448 


-2.806210 X 10-" 


2.177780 


rev 


9.507767 X 10" 


-8.064315 


-4.574355 x 10" 


-3.572965 x 10" 


3.975448 


-2.806210 X 10"" 


3.677780 



for other proton-rich Ni isotopes. The new rate is lower 
by a factor of about 0.37 than the previously used stan- 
dard rate, but the change has no significant influence on 
the I'p process calculations. A generalization to all pro- 
ton captures in the vp process (or other nucleosynthesis 
processes), however, cannot be made as the deviations 
between experiment and theory are specific to the nu- 
clides investigated here. 



0822648 (Joint Institute for Nuclear Astrophysics), the 
EUROCORES EuroGENESIS research program, the 
ENSAR/THEXO European FP7 program, the Hungar- 
ian Academy of Science, and by the DOE Topical Col- 
laboration "Neutrinos and Nucleosynthesis in Hot and 
Dense Matter" under contract DE-FG02-10ER41677. 



ACKNOWLEDGEMENTS 



This work was supported by the National Science 
Foundation under Grants No. PHY-1102511, PHY- 



[18] 



M. Arnould and S. Goriely, Phys. Rep. 384, 1 (2003). 
T. Rauscher, Proceedings of Science PoS(NIC XI)059 

(2011) . [19 
S. E. Woosley and W. M. Howard, Astrophys. J. Suppl. 
36, 285 (1978). [20 
M. Rayet, M. Arnould, M. Hashimoto, N. Prantzos, and 

K. Nomoto, Astron. Astrophys. 298, 517 (1975). [21 
T. Rauscher, A. Heger, R. D. Hoffman, and S. E. 
Woosley, Astrophys. J. 576, 323 (2002). [22 
S. E. Woosley and A. Heger, Phys. Rep. 442, 269 (2007). 
W. M. Howard, B. S. Meyer, and S. E. Woosley, Astro- [23 
phys. J. Lett. 373, L5 (1991). 

W. M. Howard and B. S. Meyer, in Proc. 2nd Int. Conf. [24 
Nuclei in the Cosmos, eds. F. Kappeler, K. Wisshak (lOP 
Pubhshing, Bristol 1993), p. 575. [25 
M. Kusakabe, N. Iwamoto, and K. Nomoto, Nucl. Phys. 
A758, 459 (2005). [26 
M. Kusakabe, N. Iwamoto, and K. Nomoto, Astrophys. [27 
J. 726, 25 (2011). 

C. Travagho, F. K. Ropke, R. Gallino, and W. Hille- [28 
brandt, Astrophys. J. 739, 93 (2011). [29^ 
E. M. Burbidge, G. R. Burbidge, W. A. Fowler, and F. [30 
Hoyle, Rev. Mod. Phys. 29, 547 (1957). 

J. Audouze and J. W. Truran, Astrophys. J. 202, 204 [31 
(1975). 

C. Frohlich, G. Martinez-Pinedo, M. Liebendorfer, F.- [32 
K. Thielemann, E. Bravo, W. R. Hix, K. Langanke, and [33 
N. T. Zinner, Phys. Rev. Lett. 96, 142502 (2006). [34 
J. Pruet, R. D. Hoffman, S. E. Woosley, H.-T. Janka, and 
R. Buras, Astrophys. J. 644, 1028 (2006). [35 
S. Wanajo, Astrophys. J. 647, 1323 (2006). [36 
T. Rauscher and C. Frohlich, AIP Conf. Proc. 1484, 73 

(2012) . ' [37 
C. Frohlich, P. Hauser, M. Liebendorfer, G. Martnez- [38 
Pinedo, E. B. F.-K. Thielemann, N. T. Zinner, W. R. 



Hix, K. Langanke, A. Mezzacappa, and K. Nomoto, As- 
trophysical Journal 637, 415 (2006). 
R. Buras, M. Rampp, H.-T. Janka, and K. Kifonidis, 
Astronomy & Astrophysics 447, 1049 (2006). 
T. A. Thompson, E. Quataert, and A. Burrows, Astro- 
physical Journal 620, 861 (2005). 

S. Wanajo, H.-T. Janka, and S. Kubono, Astrophys. J. 
729, 46 (2011). 

C. Frohlich and T. Rauscher, AIP Conf. Proc. 1484, 232 
(2012). 

C. Frohlich, T. Rauscher, X. Tang, and J. W. Truran, 
Phys. Rev. C, submitted (2012). 

"Experimental nuclear reaction data (exfor)," 
www.nndc.bnl.gov/exfor/exforOO.html (2012). 
T. Rauscher and F.-K. Thielemann, At. Data Nucl. Data 
Tables 79, 47 (2001). 

T. Rauscher, code SMARAGD, version 0.9.0s(2012). 
A. Simon, S. J. Quinn, A. Spyrou, et at, Nucl. Instr. 
Meth. A703, 16 (2013). 

T. Rauscher, Astrophys. J. Suppl. 201, 26 (2012). 
W. Hauser and H. Feschbach, Phys. Rev. 87, 366 (1952). 
J. A. Holmes, S. E. Woosley, W. A. Fowler, and B. A. 
Zimmerman, At. Data Nucl. Data Tables 18, 305 (1976). 
T. Rauscher, F.-K. Thielemann, and K.-L. Kratz, Phys. 
Rev. C 56, 1613 (1997). 

T. Rauscher, Int. J. Mod. Phys. E 20, 1071 (2011). 
W. A. Fowler, Quart. J. Roy. Astron. Soc. 15, 82 (1974). 
T. Rauscher, P. Mohr, I. Dillmann, and R. Plag, Astro- 
phys. J. 738, 143 (2011). 

T. Rauscher, Astrophys. J. Lett. 432, 326 (2012). 

T. Rauscher and F.-K. Thielemann, At. Data Nucl. Data 

Tables 75, 1 (2000). 

A. Lejeune, Phys. Rev. C 21, 1107 (1980). 

L. van Wormer, J. Goerres, C. Iliadis, M. Wiescher, and 

F.-K. Thielemann, Astrophys. J. 432, 326 (1994). 



12 



[39] J. P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. 

C 16, 80 (1977). 
[40] E. Bauge, J. P. Delaroche, and M. Girod, Phys. Rev. C 

58, 1118 (1998). 
[41] E. Bauge, J. P. Delaroche, and M. Girod, Phys. Rev. C 

63, 024607 (2001). 
[42] G. Gyiirky, G. G. Kiss, Z. Elekes, Z. Fiilop, E. Somorjai, 

and T. Rauscher, J. Phys. G 34, 817 (2007). 
[43] G. G. Kiss, T. Rauscher, G. Gyiirky, A. Simon, Z. Fiilop, 

and E. Somorjai, Phys. Rev. Lett. 101, 191101 (2008). 
[44] T. Rauscher, G. G. Kiss, G. Gyiirky, A. Simon, Z. Fiilop, 

and E. Somorjai, Phys. Rev. C 80, 035801 (2009). 
[45] T. Rauscher, Phys. Rev. C 78, 032801(R) (2008). 
[46] T. Rauscher, Phys. Rev. C 81, 045807 (2010). 
[47] M. E. Scvior, L. W. Mitchell, M. R. Anderson, C. W. 

TingwcU, and D. G. Sargood, Austral. J. Phys. 36, 463 

(1983). 

[48] C. I. W. Tingwell, V. Y. Hansper, S. G. Tims, A. F. 
Scott, A. J. Morton, and D. G. Sargood, Nucl. Phys. A 



496, 127 (1988). 
[49] G. A. Krivonosov, O. I. Ekhichcv, B. A. Ncmashkalo, 

V. E. Storizhko, and V. K. Chirt, Izv. Rossiiskoi 

Akademii Nauk, Ser. Fiz. 41, 2196 (1977). 
[50] C. I. W. Tingwell, V. Y. Hansper, S. G. Tims, A. F. Scott, 

and D. G. Sargood, Nucl. Phys. A 480, 162 (1988). 
[51] "National nuclear data center, nudat 2 database," 

http:/ /www. nndc.bnl.gov/nudat2 (2010). 
[52] "Evaluated nuclear structure datafile (ensdf)," 

http://www.nndc.bnl.gov/ensdf (2010). 
[53] G. G. Kiss, T. Sziics, Z. Torok, Z. Korkulu, G. Gyiirky, 

Z. Halasz, Z. Fiilop, E. Somorjai, and T. Rauscher, Phys. 

Rev. C 86, 035801 (2012). 
[54] C. W. Tingwell, L. W. Mitchell, M. E. Sevior, and D. G. 

Sargood, Nucl. Phys. A 439, 371 (1985). 
[55] C. W. Cheng and J. D. King, Canadian J. Phys. 58, 1677 

(1980). 



